COMBINATORIAL APPROACHES TO HOPF BIFURCATIONS IN SYSTEMS OF 

INTERACTING ELEMENTS 



DAVID ANGELI\ MURAD BANAJI^, AND CASIAN PANTEA^ 

Abstract. We describe combinatorial approaches to the question of whether families of real matrices admit 
pairs of nonreal eigenvalues passing through the imaginary axis. When the matrices arise as Jacobian matri- 
ces in the study of dynamical systems, these approaches provide necessary conditions for Hopf bifurcations 
to occur in parameterised families of such systems. The techniques depend on the spectral properties of ad- 
ditive compound matrices: in particular, we associate with a product of matrices a signed, labelled digraph 
termed a DSR'^' graph, which encodes information about the second additive compound of this product. A 
condition on the cycle structure of this digraph is shown to rule out the possibility of nonreal eigenvalues 
with positive real part. The techniques developed are naturally applicable to systems of interacting elements 
termed "interaction networks" , of which networks of chemical reactions are a special case. 



1. Introduction 

The material in this paper can be motivated both by abstract questions in hnear algebra about the spectra 
of sets of matrices, and by the study of asymptotic behaviour in dynamical systems. The connection is quite 
natural: given a sufficiently smooth dynamical system, the structure or spectra of the Jacobian matrices 
associated with the system may determine certain behaviours, for example the possibility of various local 
bifurcations, and more generally the possibility of multiple steady states, oscillations or chaos. Given a set 
X and a map J : X ^ R"^", define 

J = {J{x) : x(^X} and Spec J = |J Spec J(a;) . 

We may ask, for example: 
Ql. "Is Og Spec J?" 

Q2. "Does Spec J'\{0} intersect the imaginary axis?" 

Q3. "Does the nonreal part of Spec J7 intersect both left and right open half-planes of C?" 
and so forth. The best-known special case is where ^7 is a "qualitative class", namely it consists of all 
matrices with some sign pattern, in which case Ql reduces to the well-studied combinatorial problem of 
characterising sign nonsingular matrices [1]. By an easy argument, given convex open Y C R", and a 
function j -.Y ^ R", sign nonsingularity of the Jacobian matrices D/ in fact implies injectivity of / [2]. 
Apart from such applications, exploration of the zero patterns of sign nonsingular matrices has led to a rich 
combinatorial theory [3j. 

If J arises as the Jacobian matrices of some family of vector fields, a negative answer to either Q2 or Q3 
implies that this family does not admit Hopf bifurcation. In this paper it will be Q3 which is of primary 
interest. 

The philosophical approach taken is quite analogous to that taken in the study of sign nonsingularity. As a 
general principle, when posing some question about a set of matrices J , one hopes to associate with elements 
of J discrete objects such as graphs which are constant or vary little over J ^ and then reduce the question 
to some finite computation on these objects. The most obvious example is again when J7 is a qualitative 
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class, which can naturaUy be associated with a signed digraph, or a signed bipartite graph. However, more 
generaUy, products of matrices which have constant sign pattern (or some of which are constant) also arise 
naturally in applications. To see why it is worth going beyond qualitative classes, we introduce the notion 
of an "interaction network" as described in [4^. 

1.1. Interaction networks. Consider a system consisting of n species with values (amounts, concentrations, 
populations, etc.) xi, . . . , a;„ G X C M. Define x — (xi, . . . , XnY G to be the state of the system. Suppose 
that amongst these species there are m interactions, each involving some nonempty subset of the species, 
and occurring at rates Vi : X" R (i = 1, . . . , m) dependent on the state x, but independent of time. Define 
the rate function t; : X" ^ by v{x) = (wi(a;), . . . , u™(x))*. Finally, for i = 1, . . . , n, let /, : M™ R 
describe the total rate of change of species i as a function of the interaction rates. The evolution of the 
system is then given by: 

(1) Xi = fi{v{x)), i = l,...,n. 
or more briefly, defining / : R'" ^ R" by /(y) = (/i(y), . . . , /„(y))*, 

(2) i = fiv{x)). 
Assume that / and v are functions, so that by the chain rule, 

(3) D{f{v{x)) = Df{v{x))Dv{x). 

For our purposes here, the important point is that the right hand side of an interaction network ([2]) is 
defined as a composition of functions, and hence its Jacobian matrices have a natural factorisation ([3]). This 
factorisation was the starting point for quite general combinatorial approaches to Ql for interaction networks 
in [4], and the treatment here relies heavily on the ideas first presented there. 

Remark 1.1. (0) is not restrictive since an arbitrary system of autonomous first order differential equations 
can be represented in this way by choosing, for example, one ofv or f to be the identity function. Thus casting 
a system as an interaction network cannot be done uniquely, and can be seen as a formalism for studying 
certain questions about the system, rather than a categorisation of the system. However the usefulness of 
this construction is most apparent when f and v are defined by natural physical constraints as we shall see 
below. 

1.2. Chemical reaction netv^rorks. O.D.E. models of systems of chemical reactions, termed "chemical 
reaction networks" or CRNs, naturally take the interaction network form: in fact ([2]) reduces to 

(4) X = Tv{x). 

where F is now a constant n x m matrix, termed the "stoichiometric matrix" of the system, and v{x) is the 
vector of reaction rates. As the state variables are chemical concentrations, the natural state space is the 
nonnegative orthant in R", i.e. x € R>q. The Jacobian matrices become 

rDv{x). 

Moreover, under mild physical assumptions, the sign pattern oi Dv may be related to that of F. For example, 
if X G intR"^, no chemical occurs on both sides of any reaction, and all chemical reactions are reversible, 
then we expect Dvji to have the same sign as — F^j. This will be discussed further below. 

1.3. A roadmap for the paper. The results will be developed as follows. In Section [5] a considerable 
volume of preliminary material is gathered for subsequent use. This includes definitions, notation and 
known results, along with a few new but relatively straightforward lemmas. Section |3] constructs an object 
termed the "DSR^^l graph", which is the central tool required for everything to follow. This is followed, 
in Section |4] by a series of claims which can be made immediately from observation of the DSR^^^ graph 
combined with the results in Section [21 and without the need for deeper theory. Section [5] contains technical 
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results on the DSR^^l graph which allow, in Section |6] some nontrivial extensions of the results of earlier 
sections. Section [6] also describes limitations of the theory developed, and suggests the way forward. 

Examples are interspersed throughout, and we return to examples to note how the analysis becomes more 
sophisticated or more rapid as the theoretical tools are developed. Similarly we sometimes prove results 
which later become corollaries of deeper results; although this adds to the length of the paper, it is intended 
to make the material more transparent. Remarks linking results here to related published work, and also 
to the study of CRNs, are interspersed throughout. The frequent reference to CRNs reflects the fact that 
although these systems are only one possible application for the theory, they were an important motivator 
for this work. 



2.1. Some background on matrices. This paper is concerned with sets of real matrices. These can be 
naturally identified with subsets of Euclidean space, thus inheriting topological properties such as openness, 
closedness and connectedness. We begin with the most fundamental of these sets: a matrix M e j^nx™ 
determines the qualitative class Q(M) C jR^x™ consisting of all matrices with the same sign pattern as M, 
i.e., X € Q(M) if and only if (M^^- > 0) ^ {Xi^ > 0); {Mij < 0) ^ {Xij < 0); and (My = 0) ^ (Xy = 0). 
The closure of Q(M) will be written Qo(M), while the closure of a more general C K"^™ -^jn be written 
A4. A square matrix M is sign nonsingular if all matrices in Q{M) are nonsingular and sign singular if 
all matrices in Q(M) are singular [I]. 

Remark 2.1. Jacobian matrices of CRNs. Returning to ^ and the subsequent discussion, the notation 
developed above allows us to abbreviate the condition "Dvji has the same sign as —Tij for all as Dv G 
Q(— r*). More generally, for systems of reactions where no chemical occurs on both sides of any reaction, 
we expect the condition 



to hold on the entire nonnegative orthant, including the boundary. Even when this does not hold (namely 
when some chemical occurs on both sides of some reaction), provided every reversible reaction is treated as 
a pair of irreversible reactions, we still expect Dv to belong to the closure of a qualitative class which can be 
inferred from the network structure. Thus quite generally, the Jacobian matrices of a CRN take the form of 
a constant matrix times a matrix belonging to some qualitative class. 

2.1.1. Submatrices, minors and terms in a determinant. Given an n x m matrix M, and nonempty a C 
{1, . . . , n}, (3 C {1, . . . , m}, define M(a|/3) to be the submatrix of M obtained by choosing rows of M from 
a and columns from /3, and (provided \a\ — |/3|) Af[ck|/3] to be the corresponding minor of M. M[a\ will 
be an abbreviation for the principal minor A/[a|Q!]. Given an n x n matrix M and some permutation a of 
{1, . . . ,n} denote by P(cr) the sign of this permutation, and by the term P{a)Il2^iMic- in detM, so 
that det M = J2a ^^(<t) • 

2.1.2. The Cauchy-Binet formula. Consider an n x k matrix A and an fc x m matrix B for arbitrary positive 
integers n, m, k. The Cauchy-Binet formula ([5] for example) tells us that for any nonempty a C {1, . . . , n}, 



2. Preliminaries 



(5) 



Dv e Qo(-r*) 



/? C {l,...,m} with \a\ = |^|: 



(6) 





7 C { 1 , . . . , m } 

|7| = |a 



This formula for the minors of a product is central to many of the matrix and graph-theoretic results 
developed and cited below. This formula also has a briefer statement in terms of compound matrices 
mentioned later. 
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2.1.3. Matrix spectra and stability. Given a square matrix M, define Spec(M) to be the multiset of eigenval- 
ues of M. Denote by C_ (resp. C_) the open (resp. closed) left half-plane of C with C+ and C+ similarly 
defined. A square matrix AI whose spectrum lies entirely in C+ (resp. C+) is positive stable (resp. posi- 
tive semistable). M is Hurwitz if all of its eigenvalues lie in C-_, which occurs iff —AI is positive stable. 
P-matrices are square matrices all of whose principal minors are positive. Po-matrices are matrices in 
the closure of the P-matrices, i.e., all of whose principal minors are nonnegative. Here a P-matrix (or Pq- 
matrix) will refer to a real matrix. Of particular importance will be the following restrictions on the spectra 
of Po-matrices and P-matrices proved in |6, : 

Lemma 2.2. A complex number A — re*^ is an eigenvalue of an n x n P^-matrix if and only if 
(7) |0 - 7r| > Ti/n. 

A is an eigenvalue of an n x n P-matrix if and only if 

— 7r| > 7r/n . 

This implies, in particular that no Pp-matrix can have a negative real eigenvalue. Obviously, if — J is a 
Pp-matrix, then J cannot have a positive real eigenvalue. We will also need the following well-known fact 
about the spectrum of a product of matrices: 

Lemma 2.3. Given A G R"^"',P G R"'^", Spec(AP)\{0} = Spec(PA)\{0}. 

2.1.4. Compound matrices. We define A'^M" as the fcth exterior power of M". In particular, A°M" = M, 
A^M" = M", and for k > 2, the elements of A'^R" are finite formal linear combinations of elements of the 
form 

Ml A U2 A . . . A Ufc 

where Ui £ R", and "A" is the wedge-product (see [2 for example). Let M denote a linear transformation 
from to M", and also its matrix representation in some basis. M determines transformations from A'^R™ 
to A'"'R" in two important ways [5]: 

(1) AfC^') : A'-^M™ ^ A'^'M" is defined by: 

M^''\ui A • • ■ A life) = (Mui) A • • • A (Muk) . 

In this case the map M'^'^' is termed the fcth exterior power of M, or the fcth multiplicative compound 
of M. Choosing some bases for A'^R™ and A'^R", the resulting (^') x (™) matrix will also be denoted 

(2) If TO = n, Af ['=1 : A'^'R" ^ A'^R" is defined by: 

k 

MI'^I (mi a • • • a Ufc) = ^ Ml a • • • (Mu,) A • • • A Mfc . 

i=l 

In this case the map Afl'^l is termed the fcth additive compound of Af. Choosing a basis for A'^'R", 
the resulting matrix is (^) x (^), and will also be denoted Afl'^l Note that (-Af)!^] = -(Afl^l) go 
we can write without ambiguity — Af^^l. 

Remark 2.4. Multiplicative and additive compound matrices can easily be shown to satisfy the relationships 
(AB)^''^ = A^'^^B^'^^ (provided A and B are of dimensions such that AB makes sense) and {A + B)^^"^ = 
^[fc] ^ (provided A and B are square and of the same dimensions) . (AP)''^-' = ^(*^)p('=) is in fact a 
restatement of the Cauchy-Binet formula in terms of multiplicative compounds. 
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2.1.5. Spectra of compound matrices. Compound matrices of a square matrix M G M"^" appear naturally 
in the study of various questions on the stability of differential equations (see [9] for example) . We will refer 
to some of this theory, but most important here are basic spectral properties of compound matrices: for 
fc = 1, . . . , n, the (2) eigenvalues of M^'^^ are precisely products of all sets of k distinct elements of Spec(M), 
while the (^) eigenvalues of Aft*^! are the sums of all sets of k distinct elements of Spec(M). Here it is second 
additive compound matrices which will be of most interest: if Spec(A/) = {Ai, . . . , An} then 

Spec(M[2l) = {A, + Xj : I <i < j <n}. 

Second additive compound matrices are sometimes presented by means of a bialternate product . This 
reference also demonstrates that second additive compounds are natural objects to consider when studying 
Hopf bifurcation. For us, the following two observations are important: 

Lemma 2.5. If M is a real square matrix and is nonsingular, then SpecAf\{0} does not intersect the 
imaginary axis. 

Proof. If Specil/\{0} intersects the imaginary axis, then M has a pair of eigenvalues ztzw where a; 7^ 0. 
These sum to zero and hence M^^l is singular. □ 

Lemma 2.6. If M is a real square matrix and M^^l is a Po-matrix, then the nonreal part of Spec M (i.e. 
SpecM\K; does not intersect If-M^^^ is a Po-matrix, then the nonreal part of Spec M does not intersect 
C+. 

Proof. If the nonreal part of SpecM intersects C_, then M has a pair of eigenvalues azLiuj where a < 0. 
These sum to 2a < and hence M^^l has a negative real eigenvalue, forbidden by Lemma [2.21 The second 
statement follows immediately since the spectrum of — Af is simply the reflection in the imaginary axis of 
that of M [2]. □ 

Remark 2.7. // Af is a PQ-matrix, then in fact Lemma \2.6\ extends immediately as follows: given any 
nonnegative diagonal matrix D. the nonreal part of Spec {M + D) does not intersect C_ . This follows because 
(i) (A/ + Z?)[^l = A/[^1+Z?[^l, (ii) IJl^l is itself a nonnegative diagonal matrix, and (Hi) the set of Pq -matrices 
is closed under the addition of nonnegative diagonal matrices. This fact is of importance in applications to 
CRNs, where diagonal terms represent arbitrary degradation or outflow of chemicals (see |lHll2j ). This will 
be commented on further below. 

We also have the following two lemmas about stability: 

Lemma 2.8. Consider a square matrix M. If both M and A/[^l are P^-matrices (resp. P-matrices), then 
M is positive semistable (resp. positive stable). 

Proof. Since Af is a Po-matrix (resp. P-matrix) by Lemma [2.21 the real part of SpecAf does not intersect 
C_ (resp. C_). Since Af'^' is a Po-matrix (resp. P-matrix), by Lemma the nonreal part of Spec A^ does 
not intersect C_ (resp. C_). Combining these two facts, Af is positive semistable (resp. positive stable) . □ 

Remark 2.9. Continuing from Remark \2. 7[ if the conditions of Lemma \2.8\ hold, namely both M and Afl^l 
are Po-matrices (resp. P-matrices), then given any nonnegative diagonal matrix D, M -\- D is positive 
semistable (resp. positive stable) 

Lemma 2.10. Consider a set Ai C R"^" which is path- connected. Suppose that for each M £ Ai, both M 
and A/[^l are nonsingular and that Ai includes a positive stable (resp. Hurwitz) matrix. Then Ai consists 
entirely of positive stable (resp. Hurwitz) matrices. 

Proof. Let M' e 7W be a positive stable (resp. Hurwitz) matrix which exists by assumption. Suppose Ai 
includes a matrix M" which fails to be positive stable (resp. Hurwitz). Since Ai is path-connected, and 
eigenvalues of a matrix depend continuously on the entries of a matrix, there must lie, on any path in Ai 



6 



connecting M' and Af", a matrix M G Ai with some eigenvalue on the imaginary axis. However, since both 
M and M[21 are nonsingular, by Lemma [2751 M has no eigenvalues on the imaginary axis. □ 

Remark 2.11. The stability criterion in Lemma \2.1U\ is closely related to that of Li and Wang in [13,, and 
indeed, it is possible to use the spectral properties of compound matrices to develop other stability criteria 
for sets of matrices which can be more useful in practice than attempting directly to check the Routh-Hurwitz 
conditions. 

2.2. Hopf bifurcations. A one-parameter family of vector fields f{x, ^) generically undergoes a Poincare- 
Andronov-Hopf bifurcation (from here on abbreviated to "Hopf bifurcation") if there exists some (xc/io) 
such that /(a;o, /io) = 0, Df{xo, /io) has a pair of imaginary eigenvalues and moreover these eigenvalues cross 
the imaginary axis as /i varies through /ip ( |14j for example). Thus vector fields can be shown to forbid Hopf 
bifurcation either if the nonzero spectrum of each Jacobian matrix avoids the imaginary axis, or the nonreal 
spectrum of each Jacobian matrix avoids C_ or C+. 

Define 5R and 3 to be the real and imaginary axes in C. As indicated earlier, we are particularly interested 
in the case where the matrices factorise. In particular, when studying Hopf bifurcation, we are led to the 
following questions: 

(1) Given A C R"x™ and B C R'"^^", is 

Spec(AB)\{0} n 3 = for all A £ A S e B ? 

(2) Given A C R"x™ and B C M™^", is 

Spec(AB) C » U C7 for aU A e B e B ? 

Since under certain assumptions, the Jacobian matrix of a CRN is of the form —AB where B E Qo{A*) (see 
Remark 1 2. II) . when analysing the behaviour of CRNs we are led to the questions: 

(1) Given A e R"^^™, is 

Spec(ylB)\{0} n 3 = for all B e Qo{A^) ? 

(2) Given A e R"><", is 

Spec(AB) C SR U C7 for ah B e Qo{A') ? 
In this paper we begin the process of developing tools which provide answers to these questions. 

Remark 2.12. From Lemma \2.6l the reader may already guess that our main tools for ruling out Hopf 
bifurcations here will involve proving that jl^l or — jl^l is a Pq -matrix for each allowed Jacobian matrix of 
some system. Continuing the theme of Remark \2.7\ if Hopf bifurcation is ruled out for a chemical system 
with Jacobian matrices J because — jl^l is a Po-matrix for all J G J, then it is in fact ruled out even if 
arbitrary degradation or outflow reactions are added. (Analogous statements follow replacing jl^l " with 
"—J" and "Hopf bifurcation" with "saddle-node bifurcation", but this is well-studied [¥lll5j and not the theme 
of this paper.) 

As indicated earlier, second additive compounds arise naturally in the study of Hopf bifurcations, a theme 
discussed thoroughly in [10]. In that work, the focus lies on devising continuation methods to detect Hopf 
bifurcations, with jl^^ occurring in algebraic systems which augment the equilibrium condition. On the 
other hand, our goal here is to present conditions which rely on combinatorial features of jl^l to forbid 
Hopf bifurcation. These conditions translate into graphical requirements involving DSR graphs, which are 
discussed below. 
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2.3. The DSR graph. A class of generalised graphs, sometimes termed SR graphs and DSR graphs, have 
become useful tools in the study of properties of interaction networks, such as multistationarity and stability 

DSR graphs for general interaction networks were constructed originally in [^. The definition presented 
here involves some minor modifications of that in [4] for ease of presentation and maximum generality. Given 
an n X TO matrix A and an m x n matrix B, the DSR graph Ga,b is defined as follows: it is a signed bipartite 
digraph on n + to vertices. Si, . . . , Sn, termed S-vertices, and . . . , termed R- vertices, with arc Rj Si 
if, and only if, Aij ^ and arc SiRj if, and only if, Bji ^ 0. The arc RjSi is said to have R-to-S orientation 
and is given the sign of Aij; the arc SiRj is said to have S-to-R orientation and is given the sign of Bji. In 
addition: 

(1) A pair of antiparallel edges of the same sign are represented as a single undirected edge, so that each 
edge either has S-to-R orientation or R-to-S orientation or both, in which case we say that the edge 
is undirected. When discussing the degree of a vertex, an undirected edge incident into the vertex 
contributes exactly 1 to the degree. 

(2) It is useful to give the (directed or undirected) edge RjSi associated with a nonzero entry Aij the 
label \Aij\. We refer to this label of an edge e as /(e). If e has only S-to-R orientation (i.e. it is 
associated only with an entry from B), then we formally define /(e) = oo. This choice of labelling 
convention springs from the fact that we can generally consider A to be constant if the system being 
studied is a CRN (see Remark l2.ip . 



Remark 2.13. What is here termed G a, b would be Ga,b<: if^ the terminology off^ . Further, the construction 
here is presented for a single pair of matrices rather than directly for sets of matrices; note however that if 
A is fixed and B varies over a qualitative class, then the DSR graph remains constant. On the other hand, 
if A varies (within a qualitative class), then varying edge labels can be replaced with the formal label oo. 

Example 2.14. Below we illustrate the DSR graph Ga,b associated with two matrices A and B. Negative 
edges ( corresponding to negative entries in the matrices ) are shown as dashed lines, while positive edges 
are bold lines, a convention which will be followed throughout this paper. To make the construction more 
transparent, an intermediate stage is shown where separate signed digraphs Ga.o Gq.s o,re constructed 
from A and from B: Ga,b can be regarded as the "superposition" of these two digraphs. Edge-labels are 
associated only with A: an edge imported only from B is labelled oo. 



Si ■* — i?2 ~ — *■ *S'2 



A = 




GAfi 



B' = 




Gq.b — 



Ri 



Si 



Ri 



S. 



■R2 



Q 3 D ^ Q 



So 



S3 




J 



Ga, 



Remark 2.15. A special case relevant to the study of CRNs is when the matrices A, B satisfy B e Q(^*). 
In this case all edges of Ga,b oltc undirected, and in fact Ga,b = Ga,a* ■ This graph can also be referred to 
as the SR graph corresponding to the matrix A. 



To present the main results about DSR graphs, some further definitions are needed. 
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2.3.1. Walks and cycles. A walk W in a digraph is an alternating sequence of vertices and edges, beginning 
and ending with a vertex, and where each edge in W is preceded by its start-point and followed by its 
end-point, each vertex (except the first) is preceded by an edge incident into it, and each vertex (except the 
last) is followed by an edge incident out of it. Here we allow a walk W to repeat both vertices and edges, 
and its length \W\ is the number of edges in W, counted with repetition. If the first and last vertex are 
the same, the walk is called closed. In what follows we may refer to walks by their sequence of edges, or 
vertices, or both and we will say that a walk is empty and denote it as if it contains no edges (an empty 
walk may include a single vertex). When used without qualification, "walk" will mean a nonempty walk. 
The sign of a walk is defined as the product of signs of the multiset of edges in the walk. The empty walk 
is formally given the sign +1. 

A closed nonempty walk which does not repeat vertices (except in the trivial sense that the first and last 
vertices are the same), is called a cycle. Clearly a cycle has no repeated edges. A DSR graph with no cycles 
will be termed acyclic: note that the walk (m, w, u) is considered a cycle in a DSR graph if and only if arcs 
uv and vu have different signs and so are not treated as a single edge. 

2.3.2. Parity of walks, e-cycles and o-cycles. As DSR graphs are bipartite objects, each closed walk W has 
even length. Any walk of even length can be given a parity P{W) defined as: 

P(W) = (-l)l^l/2sign(iy). 

If P{W) — 1, then W is termed even, and if P{W) = —1, then W is termed odd. A cycle C is termed an 
e-cycle if P{C) ~ 1, and an o-cycle if P(C) = — 1. 

2.3.3. Closed s-walks and s-cycles. A closed walk W = (ei, 62, ... , e2r) is an s-walk if each edge einW has 
an associated label ^(e) ^00 and 

r r 

l[l{e2^-l)=l[l{e2^)■ 

i=l i=l 

If is a cycle, it is termed an s-cycle. 

2.3.4. Intersections of cycles. A cycle in a DSR graph which includes some edge with only S-to-R orientation 
or only R-to-S orientation has a unique orientation; otherwise the cycle has two possible orientations. Two 
oriented cycles in a DSR graph are compatibly oriented, if each induces the same orientation on every edge 
in their intersection. Two cycles (possibly unoriented) are compatibly oriented if there is an orientation for 
each so that this requirement is fulfilled. In a DSR graph, two cycles have odd intersection if they are 
compatibly oriented and each component of their intersection has odd length. Note that odd intersection 
was termed "S-to-R intersection" in [TCI H] . 

Remark 2.16. A necessary condition for the existence of two cycles, say C and D, with odd intersection 
in a DSR graph is the existence of at least one S-vertex of degree 3 or more, and of at least one R-vertex of 
degree 3 or more. This follows since the vertices at the ends of a component of the intersection of C and D 
of odd length must have degree at least 3, and one of these must he an S-vertex and one an R-vertex. 

2.3.5. Hatchings. A matching in a DSR graph is a set of edges without common vertices. An S-to-R matching 
is a matching all of whose edges have S-to-R orientation, and an R-to-S matching is similarly defined. A 
matching which covers all vertices is perfect. 

2.4. Some useful results on matrix sets and DSR graphs. We will refer to a DSR graph G as odd if 
G has no e-cycles. We will refer to it as odd* if all e-cycles in G are s-cycles, and no two e-cycles have odd 
intersection. It will be called steady if all cycles in G are s-cycles. We summarise several important results 
connecting DSR graphs and Pi^-maXiices in a single theorem. 

Theorem 2.17. Define the following conditions on matrices A G ][j"x™ and B G M™^" and associated 
graphs: 
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CI. Ga.b is odd. 
C2. Gam is odct . 

C3. A'B' is a Pq -matrix for all A' e QoiA), B' e Qq{B). 
C4. AB' is a Pq -matrix for all B' e Qo{B). 
C5. Ga^a^ is odd. 
C6. GA,At is odd* . 

C7. Every square submatrix of A is either sign nonsingular or sign singular. 
C8. Every square submatrix of A is either sign nonsingular or singular. 
C9. A'C is a P^-matrix for all A' e Qo(^), C £ Qo{A*). 
CIO. AC is a Po-matrixfor all G G Qo(^*)- 

The following implications hold: 





CI 




C3 
















C2 




C4. 




C5 




C7 




C9 












C6 




C8 


<^ 


CIO 



Proo/. The implications CI =^ C2, C3 C4, C5 ^ C6, C7 =^ C8 and C9 => CIO follow by definition, 
while the other implications are results in [151 [12111], or are immediate consequences of these results. Note 
that C5 <^ C9 and C6 ^ CIO are just immediate corollaries of CI <^ C3 and C2 => C4 respectively. □ 

Remark 2.18. Although some of the proofs of results in Theorem \2.17\ are lengthy, the spirit is worth 
outlining. Consider matrices A, B and DSR graph Ga.b- Since the results are about the product AB, they 
begin with the Cauchy-Binet formula which provides an expression for each minor of AB as a sum of 
products of minors of A and B. The transition from matrices to DSR graphs is via the association of terms 
in subdeterminants with matchings in DSR graphs. Given appropriate sets a, (3 a term in A[a\l3] corresponds 
to an R-to-S matching in Ga.b which is perfect on the vertices {Si}i£a U {-Rjjjg^, while a term in B[l3\a\ 
corresponds to an S-to-R matching which is perfect on the same vertex set. The union of these matchings 
is a subgraph of Ga.b consisting of (disjoint) cycles and isolated, undirected, edges. Conclusions about the 
signs and values of terms in the minors follow from examination of this subgraph. 

Alongside Theorem 12 .171 we will need a few further results specifically about steady DSR graphs, which 
were not considered in the references above, but become important in the discussion here. 

Lemma 2.19. Let G be a steady DSR graph. Then any closed walk on G is an s-walk. 

Proof. Let be a closed walk on C. Certainly each edge of W which occurs in a cycle of C has finite label, 
since all cycles are s-cycles. Further, each edge of W not occurring in a cycle of C must be undirected, and 
consequently have a finite label: traversing arc uv and returning subsequently to u otherwise implies the 
existence of a cycle which includes uv. Thus all edges in W have finite labels. We show that 

x{w)^ I n '(^)| ( n ^(^)) =1- 

\ee{S-to-R edges of W} / \ee {R-to-S edges of W} / 

If no vertex is repeated in W , then W is either a cycle, in which case A(W^) = 1 by hypothesis, or a 
single edge traversed back and forward, in which case again \{W) = 1. If, on the other hand, W contains 
repeated vertices, then the shortest closed subwalk of W, denoted U, is either a cycle or a single repeated 
edge. Denote by W\U the closed subwalk obtained by removing U from W, namely if the vertex-list of W is 
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(wi, 102, ■ . ■ — wi) and that of U is {wj, . . . ,Wk = Wj), then for W\U we have (wi, . . . ,Wj, Wk+i, . . . ,Wr = 
wi). Clearly 

A(W^) = \{U)X{W\U) = \{W\U). 

If the walk W\U does not repeat vertices, we are done. Otherwise, we may repeat this procedure until the 
remaining walk has no repeated vertices, and the lemma is proved. □ 

Since an acyclic DSR graph vacuously fulfils the conditions of Lemma l2.19| an immediate corollary is: 
Corollary 2.20. Let G he an acyclic DSR graph. Any closed walk on G is an s-walk. 

Lemma 2.21. Consider some M' € M"^™, and some M C Q(Af') such that for each M € M, Gm,m* is 
steady. Given any A,B*gA4, AB is a Po-matrix. 

Proof. If the premise of the theorem implies that AB is a Pg-niatrix for each A,B* G Al, then the result 
follows for A, i3* G A4 since the set of Po-matrices is closed in R"^"-. So assume that A, B* G Ai. Given any 
nonempty a C {1, . . . , n}, by the Cauchy-Binet formula: 



a 



(8) iAB)[a]^ A[ah]B[j\ 

7C{ 1 m} 

|7| = |a| 

Fix 7, define M — A(a\j), and k — \a\. Let M^g.) and be two nonzero terms in detM = ^[q!|7] 

corresponding to permutations a =^ t of {I, . . . ,\a\}. Note that 

k k 
i=l 1=1 

Mi-g.-) and Af(^-) each correspond to a matching in Ga.a*: and the union of these two matchings defines a set of 
disjoint cycles corresponding to nontrivial cycles in the permutation Toa~^ and isolated edges corresponding 
to trivial cycles in r o cr~^ (see Remark l2.18p . Since all cycles in Ga,a* are s-cycles, we get = |A'/(t)I 

(see the proof of Lemma 7 in [12j for example). Thus ^[a|7] = r{a,"f)\M(^„j\ where r(a,7) is an integer 
independent of the values of nonzero entries in A. Since -B* belongs to the same qualitative class as A and 
all cycles in Gst^B are s-cycles, we get -B[7|a] = r(a, 7)|-/V(ct) | where N = B{-f\a). Thus 

A[a\j]B[j\a] = [r{a,j)f |M(,)| |7V(,)| > 0, 

and the result now follows from (HI). □ 



3. Compound matrices and the DSR'^I graph 

The main purpose of this section is as follows: given a square matrix product AB, we wish to be able to 
write (Ai?)[^l = AB where A depends only on A, and B depends only on B. This will allow us to construct 
an object, the DSRI^' graph of AB, whose structure encodes information about (AB)^'^^ allowing the use of 
Theorem 1 2 . 1 71 and subsequent results. To motivate the construction we first present two examples: 

Example 3.1. Consider the matrices 

/110\ / a b \ / a + b c + d f \ 

(9) A = i 1 1 1 , B*= c d e , J = AB=\ a + b c + d + e f + g . 

V 1 1 / \0 f g J \ b d + e f + g J 

Assume that a,b,c,d,e, f,g > 0. Here, as in all examples to follow, we let Ci denote the ith standard basis 
vector in K" and choose the natural basis for A^R'^ consisting of the wedge-products of pairs of vectors Ci, 
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arranged in lexicographic order. So for A^M"^, using the basis (ei A 62, ei A 63, 62 A 63), we can compute: 

/ a+b+c+d+e f+g -f \ 

(10) Jl^l = d + e a + b + f + g c + d 

\ -b a+b c+d+e+f+gj 

By computing its 7 principal minors, we can check directly that J^^^ is a Pa-matrix, although J itself is not 
a Po-matrix. We will see later that the conclusion that j'^' is a P^-matrix can he reached without symbolic 
computation, and this example is a special case of a much more general result. 

Example 3.2. Consider the matrices A, B defined by: 



A = 



( " 


b 








a' 


b' 


' \ 


— c 


d 









-c' 


d' 








e 


/ 









e' 


/' 


V 





9 ) 




[ 








9' 1 



where a, b, c, d, e, /, g, a' , b' , c' , d' , e' , /', g' are all positive real numbers. AB is a 4 x 4 matrix, so (AB)^^^ is a 
6x6 matrix, (i) We can check that AB is a Po-matrix. (ii) With the help of a symbolic algebra package such 
as MAXIMA [20j . we can compute the 63 principal minors of {AB)^'^\ confirm that these are all nonnegative, 
and conclude that (AB)^^^ is a Po-matrix. By Lemma \2. 8\ (i) and (ii) imply that AB is positive semistable. 
It is again natural to ask whether this conclusion can be reached simply, via Theorem \2. 17\ and without the 
use of symbolic algebra. We will return to this example and show that this is indeed the case. 

3.1. Identifying A^M" with antisymmetric matrices. We will eventually write: 

„ 

where the niatrix L depends on A and L depends on B. The linear algebraic justification for the factori- 
sation begins by identifying A^K" with the set of n x n antisymmetric matrices. 
Define the matrix Xij by 

Xij = CiCj — CjCj^ . 

Note that {XijY — Xji = —Xij. Denote by asym(n) the subspace of M"^" spanned by {Xij}, i.e. the 
space of all real n x n antisymmetric matrices. Given any C € R"^" define L*-^ : asym(n) — )■ asym(n) via 
L^{X) = CX + XC*. Identifying X^j e asym(n) with A Cj G A^M", L'^ can be identified with the second 
additive compound C^^'. To be precise, define l : A^M" — >■ asym(n) via i(ei A Cj) = Xij. That this action 
on basis vectors extends naturally to give a well-defined linear bijection can easily be checked. We can then 
confirm that given any square matrix C, the following diagram commutes: 

C[21 

A^M" ^A^M" 

L L 

asym(n) > asym(n) 
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The calculation is straightforward: 

(6oC[2l)(eiAe,) 



liCci A Cj + Cj A Cej) 



= i 



Cfejefe A + t A ^ Cfejej.^ 

k k 

= CkiXkj + CkjXik 

k k 

= X] ~ ^i^k) + Ckjieicl - eke*)] 

k 

= Ceie] - ejelC* + eie]C* - Ceje* 

= CXij + XijC' = L^Xij = {L^ o i){ei A e,) . 

3.2. Factorising the second additive compound. The construction is first carried out for rank 1 matrices 

and extended to arbitrary matrices. 

Given s € M", we can define the linear map L : ffi." — >■ asym(n) via its action on coordinate vectors: 

2—1 n 

L^i = CjS* - se- = ^ eiSiej - ^ s;e;e* = ^ siXu = siXu + ^ siXu . 

I I I 1=1 

Noting that ii i < I and j < k, then e^jXuek = SijSki (where 6ij is the Kronecker delta), a matrix 
representation of L" is obtained: 



■'{j,k),i 



n 

^ sie*Xiiek - ^ sie^Xuek 



(i < k) 



l=i+l 



(11) 



n 


i-i 




r sfc 


= ^ si6ijSki 




k^jl = 




l=i+l 


1=1 




I 


{si,S2,S3,S4,y, then 










1 2 


3 


4 


(1,2) 


/ S2 -.S'l 





\ 


(1,3) 


S3 


-Si 





r _ (1,4) 


S4 





-si 


^ " (2,3) 


S3 


-S2 





(2,4) 


S4 





-S2 


(3,4) 





S4 


-S3/ 



{^ = J) 
{i = k) 
otherwise 



Row and column labels have been shown explicitly to highlight the fact that each column is associat(;d with 
an entry in s, while each row is associated with a pair of entries. If f is a i x n matrix (a row vector) then 
L" : asym(n) K" is similarly defined via its action on basis elements of asym(n): 



V^Xjk = XjkV* = {ejel - ekej)v* = ejVk - ekVj . 



13 



With j < k, we obtain the coordinate representation: 



Vk 



[IfXjk]^ = e\{e-jVk - ekVj) = Si-jVk - SikVj 



{z = k), 



otherwise. 



Thus ^ (L" )*. 

Now consider some rank 1 matrix C — sv where s, E M", and the associated operator L*^ : asym(rt) 
asym(n). Applying the definitions: 

L^X = svX + Xv's' = {Xv*)s* - s{Xvy = (L^X)s* - siJ/XY = L . 

In other words, — L'^L^. Thus representations of L'' and L" above give us a natural representation and 
factorisation of L*". 

More generally, consider any A G M"^™ and B e M™^". Defining Ak to be the A;th column of A and B'' 
to be the fcth row of B we have AB = J2k ^kB'^ . We get 

L^^X = Y^{AkB''X + X{B''Y{Akf) 

k 

= Y.^{X{B%){AkY - Au{X{B%f) 

k 



Choosing a basis for A M" we define for any n x m matrix A the matrices 



and Jj^ 



A" 



Note that 



(12) 



Given A e 



and B e 



we get: 



L 



Observe that L has dimensions (2) x mn and has dimensions mn x (2) • 

Remark 3.3. When A and B are matrices with symbolic entries, one approach to computation of any minor 
of (y4_B)[^l is by application of the Cauchy-Binet formula to I112\) . For example assuming that A is n x m 
and B is m X n, and defining a = {!,..., (2) }• 

(13) det((AB)[2i)= ^^m'^ma]. 

/3C{l,...,nm} 

l/31 = (2) 

This observation can be used in conjunction with Lemma WT^ to rule out the possibility of pairs of imaginary 
eigenvalues of AB. 

3.3. The DSR[^1 graph. Given a (square) matrix product AB we define the DSRt^l graph of the product 
AB 

to be the DSR graph of the product L L . Examination of this graph allows us to make statements about 
(AB)[^1 using, for example. Theorem 12. 171 
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3.3.1. Terminology and notation. The S-vertices of Cj^^g are the (2) unordered pairs {i,j) G {l,...,n} x 
{1, . . . , n}, i ^ j; the R-vertices of G]^ g are the mn pairs (fc, /) G {1, . . . , m} x {1, . . . , n}. For notational 
convenience we denote S-vertices {i,j) — (j, i) by ij, or ji; R-vertices {k, I) by fc'; and an edge between them 
of unknown direction and sign by {ij, fc'). From the definitions, a necessary condition for the edge {ij, fc') to 
exist is I S {i,j}- 

Remark 3.4. Continuing from Remark \2.15\ it is clear that given a matrix A and some B G Q(A*), 

gJ^'^ — cj^'^f. It is also not hard to see that if is odd*, then so is G^^ g where B £ Qo{A*) (G^^g 

[2] 

is just a subgraph of Gj^j^t)- On the other hand, if A has more than one row, and there exist i,j such that 

[21 

AijBji < 0, then G]^ ^ automatically fails Condition C2. This follows because given any k ^ i there exist 

two edges {ik,j^) in G^'^g, oppositely directed, of opposite sign, and one having label 00; so G^'^ g has an 
e-cycle of length two which is not an s-cycle. In this paper, all the examples we present will involve matrix 
pairs A, B satisfying AijBji > 0, and most often we will in fact assume B G Qq{A*). 

Example 3.5. If A = (1, 1, 1, 1)*, we can easily compute G||'^t to be: 




All edge-labels are I, so have been omitted. More generally it is clear that if v £ R" has no zero entries, 
then (disregarding edge-signs and labels) G\j [^t is a subdivision of the complete graph Kn on the R-vertices 

4. Results illustrating the use of the DSR'^' graph 

Here we develop and apply some results using DSR^^i graphs before the development of further theory. 
Some of the examples presented will later be shown to fall into families amenable to analysis, after the 
development of theory relating DSR and DSR^^' graphs. 

We start with the following rather basic result about the second additive compounds of certain rank 1 
matrices. The result can easily be proved directly, but demonstrates application of the DSR^^l graph. A 
far-reaching generalisation of this result will be presented as Theorem 16.31 

Lemma 4.1. If s,v G M", v G Qo{s) and J = sw*, then jl^l is a Po-matrix. 

Proof. If the result is true for v G Q(s), then it follows for v G Qo{s) since the set of Pg-matrices is closed. 
So we assume that v G Q{s). By Remark 13.41 G\, — G\, \,f Note that from the definition of L , row {j, k) 
{j < k) has at most two nonzero entries, the entry —Sk in column j and the entry Sj in column k. Thus the 
edges incident with S-vertex jk have labels Sj and Sk. 

i. No two cycles have odd intersection. Since no row of has more than two nonzero entries, each S- 
vertex in Gf^t has maximum degree 2, as illustrated for example by the DSR^^l graph of ExamDle l3.5l 
This rules out odd intersections (see Remark l2.16|) . 
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ii. All cycles are s-cycles. Define X — {l,...,n}, X2 = {(1, 2), (1, 3), . . . , (n — 1,"-)} (the set of all 

unordered pairs from I). Let ii,...,im € X he distinct, and ai, . . . ,am S X2 be distinct. An 
[21 

arbitrary cycle in G|, [,t must be of the form 

where ik,ik+i G ctk, and since ik,ik+i are distinct this implies that ak = {ik,ik+i)- (Here we define 
hn+i = *i-) So the sequence of labels of edges in C is 

and C is clearly an s-cycle. 
(i) and (ii) together imply that G^^\t is odd*, and hence by Theorem 12.171 that (sw*)[^l is a Po-matrix for all 

v£Qis). ' a 

Example 13.11 continued. We now return to the matrices in We will see that the DSBP'^ graph makes 
our earlier conclusions immediate, hut shows moreover that the example provides just a special case of a 
general result, Proposition \4- 1\ below. Using ill}) , we compute 

_ /l-l 01-1 010 o\ 
L'^ = -1 1 -1 1 . 
\ 0-10 1-10 1-1/ 

The DSR graphs Ga,a* fl^ft) ^.n-d Gj^\t (fight) are depicted below: 




All edge-labels are 1 and have been omitted. Consequently, all cycles in both G a.a* o-^d G are s-cycles. 
Ga,a^ includes the e-cycles (Si, 5*2, -R2, "Si) and {82, R2, S3, R3, S2) which have odd intersection and so 

fails to be odd* . R is also easy to confirm that for some choice of B ^ Q(^*), AB fails to be a P^-matrix. On 

[2] [2] 
the other hand, since all R-vertices in G^ j^t have degree less than or equal to 2, G^^ cannot have a pair 

of cycles with odd intersection, and so is odd* (see Remark \2.1b]) . By Theorem \2.17[ (Ai?)[^l is a P^-matrix 
for all B G Q{){A^), and so, by Lemma \2.6\ the nonreal eigenvalues of AB lie in C+. 

Remark 4.2. In fact, if B g Q{A*), and the matrices J = AB arise as the Jacobian matrices of a differential 
system, the form of J^^l in hlO\) has strong consequences. In particular, suppose Y 'ZK^ is open, / : y — > M"^ 
is G^ , and consider the system x ~ f{x) on Y . Let X C Y be some compact, forward invariant, set, and 
assume that for each x G X, —Dfix) is of the form in Then, by observation of UU\). Dfix)^"^^ is 

irreducibly diagonally dominant (in its columns) with negative diagonal entries, and it is possible to find a 
(single) positive diagonal matrix E such that E Df{x)^^^ E^-^ is strictly diagonally dominant for all x ^ X. 
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Thus there exists a "logarithmic norm" IJ-e{') and p < such that for all x G X, /i£;(Z?/(x)[^l) < p [21 , 
and applying the theory developed by Li and Muldowney ( \22[ 123) for example), the system x — f{x) cannot 
have invariant closed curves in X : periodic orbits whether stable or unstable are ruled out, as are homoclinic 
trajectories, heteroclinic cycles, etc. In fact every nonwandering point of the system must be an equilibrium 

m- 

Example 13.11 suggests the following much more general result: 

Proposition 4.1. (i) Consider any 3 x m (0,1,-1) matrix A and any B G Q{A*-). Then [AB)^'^^ is a 
Pg-matrix. (ii) Consider any m x 3 (0, 1, —1) matrix A and any B G Q(A*). Then the nonreal spectrum of 
(A_B)[^1 does not intersect C_. 

Proof, (i) All edge-labels in gJ^'^ are 1 and so all cycles are s-cycles. On the other hand, all R-vertices 

[21 

in Gj^g have degree less than or equal to two, and so no odd intersections between cycles are possible. 

Consequently, G^^g is odd*, and so (A_B)[^1 is a Pg-matrix. (ii) Trivially the spectrum of AB is equal to 
that of B^A*, and in turn the nonzero spectrum of is equal to that of A^B* (Lemma 12. 3p . From (i), 

(A*_B*)[^1 is a Po-matrix, and consequently the nonreal spectrum of A^B^ cannot intersect C_. Thus the 
same holds for AB. □ 

Further generalisations of Proposition 14. II will be proved as Theorem 16 . 1 1 later . 

Remark 4.3. Proposition \4. 1\ has the following interpretation in terms of the dynamics of CRMs. Consider 
a system of reactions where all entries in the stoichiometric matrix have magnitude 1 (a fairly common 
situation) and where the kinetics satisfies the condition in (0). Suppose further there are no more than three 
substrates, or no more than three reactions. Then the system is incapable of Hopf bifurcations. 

The next result which follows with little effort, involves (0, 1, —1) matrices with no more than two entries 
in each column (or row). 

Proposition 4.2. (i) Consider any n x m (0,1,-1) matrix A with no more than two nonzero entries in 
each column and any B G Qo(^*)- Then (AB)^^^ is a P^-matrix and so the nonreal spectrum of AB does 
not intersect C_. (ii) If A has no more than two nonzero entries in each row then the nonreal spectrum of 
AB does not intersect C„ . 

Proof. In both cases, as all edge-labels in G^"^ ^ are 1, all cycles are s-cycles. (i) All that needs to be confirmed 

[21 

is that all R-vertices in Gyi\g have degree less than or equal to two, and so no odd intersections between 
cycles are possible. This follows by noting that the edge (ij, fc-') exists if and only if Aik ^ 0, which (for 
fixed fc, j) can occur for a maximum of two indices i since no column of A contains more than two nonzero 

[21 . 

entries. Consequently, each vertex has degree < 2, and Gj^g is odd*. The remaining conclusions follow 
from Theorem 12 . 1 71 and Lemma l2.6l (ii) If A has no more than two nonzero entries in each row, then by (i), 
(A*i3*)[^l is a Po-matrix and consequently the nonreal spectrum of does not intersect C_. Since the 

nonzero spectrum of AB is the same as that of A'^B*, the same holds for AB. □ 

Remark 4.4. Proposition \4 . 2\ will be generalised as Theorem \6.'A where the requirement that A is a (0, 1, —1) 
matrix will be weakened to the requirement that all cycles in Ga,a* o-f^ s-cycles. CRNs where the degree of 
each S-vertex is less than or equal to two have been discussed previously in |17) . 

Example 13.21 continued. Returning to this example we have the matrix A and DSR graph Ga^a* ■ 
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A = 



f a b \ 

-c d 

e / 

\ 9 / 



Ga,a* = Ri 



Si 



d 
\ 
i?2 



/ g 
S3 -R3 S4 



Ga.a* is odd by inspection, and, so AB is a Pq- matrix for all B G Qq{A^). Following removal of vertices and 

[21 

edges which do not participate in any cycles, we get the reduced representation of G^^t ; 



/ 



23 



which can he computed to he odd*". Consequently [AB)^"^^ is a Pq -matrix for all B g Qo{A*^)- This means 
that matrices of the form AB cannot have a pair of nonreal eigenvalues in C_ . 

Example 4.5. In this example we illustrate the use of the DSR^'^^ graph to make claims about the second 
additive compound of a single square matrix ( rather than a product of two matrices). Let a, b, c, d, e, f,g, h, j > 
and consider the matrix and DSR graph: 



Sa 



I a 6 \ 
—cdgQ 
e /i 
V 3 f J 



Ga,. 



7 
Si 

\ 



Ri 



i?4 



^3 
J. 



R2 



i?3 



Here I is the 4x4 identity matrix which we note lies in Qo{A). Ga,i has exactly two cycles, one of which 
is odd and one of which is even. Following removal of vertices and edges which do not participate in any 
cycles, we get the reduced representation of G^^j: 



23 




[21 

Edge labels have been omitted for the reason thai G]^ j is odd, namely there are no e-cycles at all. Consequently 
^[2] = (^J)[2] is a Po-matrix, and matrices of the form A cannot have a pair of nonreal eigenvalues in C_. 
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Remark 4.6. As Ga,i has an e-cycle, A is not amenable to analysis using the results of [2] or [24J . The 
sufficient conditions for the absence of stable oscillation in dynamical systems developed in these papers are 
neither necessary nor sufficient for the preclusion of Hopf bifurcation. 

5. Relationships between the DSR and DSRI^' graphs 

Throughout this section we consider two matrices A e M"^™ and B G K™^". In order to understand 
better the DSR^'^l gi'aph, it is convenient to develop some theory on the relationship between structures in 

[2] 

Ga.b and The results presented below begin this process. Applications of this theory are presented 

in the final section. 

5.1. Removing leaves. Our first observation is that a certain simple operation on DSR graphs - the removal 
of pendant R- vertices and their incident edges does not alter the cycle structure of the corresponding DSR^^l 
graph. This result does not extend to the removal of pendant S- vertices. 

Lemma 5.1. Consider an n x m matrix A (m >2) whose kth column has a single nonzero entry Akg,k, so 
that vertex Rj. in Ga.a* has degree 1. Let A be the submatrix of A where the kth column has been removed. 
Then is odd* (resp. odd, resp. steady) if and only ifG^^^t is odd* (resp. odd, resp. steady). 

[21 [21 

Proof. In one direction the result is trivial: if Gj^j^t is odd* (resp. odd, resp. steady), then so is G^^^ as 

[21 

it is a subgraph of G^^^^t. 

[21 [21 

In the other direction, we show that G^^^ is an induced subgraph of Gj( obtained by removing only 

[21 [21 

leaves of G^^t, and consequently has the same cycle structure. Note that G^^^ is the induced subgraph 

[21 . . . ■ . . . ' 

of G obtained by removing all R- vertices of the form k^ and their incident edges. However, an edge in 

[21 

G j^t of the form [ij, k^) exists if and only if Aik ^ which occurs iff i = feg- Thus for each j, column k^ 
of L has a single nonzero entry (in row koj), and the corresponding vertex in G||'^t is a leaf which can be 

[21 

removed without affecting any cycles of G^^^t- Q 

5.2. Projections of DSRl^l cycles. It is useful to distinguish between two types of edges in a DSRI^' 
graph, (ij,fc"""^*'J>) and (ij, fc^'^'^'f^'^^). Recah from Section [3] that if i < j then 

Ly,fci = Ajk , Lfi ,y — Bkj, 

whereas 

[21 

This fact, in conjunction with the definition of Ga,b, implies that for any edge {ij, fc*) in G)^^ there exists 
an edge (Sj, Rk) in Ga.b having the same direction, the same label and the same sign; and similarly, for any 

[21 

edge (ij, k^) in G^^ there exists an edge {Si, Rk) in Ga,b having the same direction and label, but opposite 

[21 

sign. We have defined a mapping of the edge-set of G)^^ into the edge-set of Ga,b- 

. . [21 [21 

Definition 5.2. Let E = E{Ga.b) and E2 ~ E{Gj^g) denote the set of edges ofGA,B andGj^^ respectively. 
Similarly let V — V{Ga.b) and V2 — V{Gj^ g). The mapping tt : E2 E defined by 

TT{{ij,k=))^{S„Rk) and 7r{{ij,k')) = {S,,Rk) 
is called the projection map of DSR^^l . 
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Remark 5.3. A few immediate properties of n are worth emphasising. 

1. TT is surjective if n > 2, but not necessarily injective if n > 3. Indeed, if {Si,Rk) G E then for any 
j G {1, . . . , n}\{i} we have 7r((ij, fc^)) = {Si, Rk)- If there is more than one choice for such j then tt is not 
injective. 

2. Following the discussion preceding Definition 15.1^1 tt acts on edge-signs as follows: 

sign(e) i/e = (ij,fc™'"{*J>) 
-sign(e) ife= {ij, fc™^'^^*'^^) 



(14) sign(7r(e)) 



[21 

3. TT does not extend naturally to vertices of G)^ g. While it may he regarded as a mapping of an R-vertex 
€ V2 to vertex Rk G V , there is no such interpretation for an S-vertex ij € V2: edges in E2 incident 

[21 

with ij may be mapped to edges in E incident with either Si or to Sj . Consequently, a walk in G^^g does 
not necessarily project to a walk in Ga.b, o,s illustrated in the picture below: 



Notice also that projection may change the sign of an edge. 

[21 

Having defined the projection map, it is natural to examine the projection of cycles in g. 

[21 

Proposition 5.1. Let G — (ei, . . . , ejy) be a cycle in G)^^, and let ij be the initial vertex of ei. Then there 
exists a partition o/{l, . . . ,N} — {ii < . . . < im} U {ji < ■ • • < Jtv-a/} such that W = (7r(ei^), . . . ,'!r{eij^,)) 
and W" — (7r(ejj ), . . . , 7^(ej„_J^J)) are walks in Ga,b, o,nd such that W' and W" have different initial vertices 
and different terminal vertices, all of which belong to {Si,Sj}. In other words, either 

i. W' and W" are closed walks, one containing Si and the other containing Sj] or 

ii. W' and W" are walks, one from Si to Sj and the other from Sj to Si . 

Proof. In this proof, when referring to an S-vertex ab in G we assume that a < b. We construct simultaneously 
the two walks, proceeding inductively along consecutive S-vertices of C, as explained below. One crucial 
feature of this construction which follows from the inductive argument is that when we arrive at the S-vertex 
ab of C, then W' and W" terminate at different vertices from the set {Sa, Sb}- 

For r = 1, . . . ,N/2, we assign the pair (7r(e2r-i), 7r(e2r)) to either W or W" as follows. At stage one, 
ij is the initial vertex of ei, and we assign Si to W and Sj to W" . At stage r, let ab denote the S-vertex 
which is the initial vertex of e2r-i, and assume, without loss of generality, that Sa G W and Sb G W". The 
subsequent pair of edges (e2r-i, e2r) must take one of the two forms 

1. {ab,k'',ab') or 2. {ab,k^,a'b) . 

In case 1 we append the walk (7r(e2r-i), 7r(e2r)) — (Sb, Rk, Sb') to W" , and in case 2 we append the walk 
(7r(e2._i),7r(e2,)) = {Sa,Rk,Sa') to W. 

Since each edge is assigned to exactly one of the two walks we have obtained the desired partition. The 
remaining properties of W and W" follow immediately from the construction. □ 

Remark 5.4. 1. From the proof of Provosition \ 5.1\ it is important to underline that a pair of consecutive 
edges of C beginning and terminating at an S-vertex, are projected to a pair of consecutive edges in either W 
or W" . In other words, in the notation of Provosition \5. 11 for each r = 1, . . . , M/2, i2r = *2r-i + 1, o,nd for 
each r — 1, . . . , {N — M)/2, j2r = J2r-i + 1- In particular, a sequence {Sa, Rk, Sa) traversing the unoriented 
edge {Sa,Rk) in- each direction consecutively cannot occur in W , W" or W' U W" : such a sequence must be 
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the projection of a sequence of edges in C of the form {ab, k^, ah) where {ab, k^) is unoriented and traversed 
hack and forward, a contradiction. 

2. In Case (i) of Provosition [XT] W and W" are closed walks in Ga.b cind it is not hard to see that 

[2] 

choosing a different initial vertex from C in Gj^ ^ leads to the same two closed walks (albeit possibly renamed, 
and traversed from different initial points) ; in Case (ii) of Provosition \5J\ W' and W" together form a closed 
walk of length N , which we will denote by W U W" . Again choosing a different initial vertex from C leads 
to the same object, traversed from a different initial point. 

[21 

It is transparent from the previous remark that whether projection of a cycle C in Gj^ g leads to Case 

(i) or Case (ii) in Proposition 15.11 is a property of C alone. A cycle C of G^^^ will be called direct if it 
projects as in Case (i) and twisted if it projects as in Case (ii). In view of Remark 15.41 the following notion 
is well-defined: 

rn] 

Definition 5.5. Let C be a cycle in Gj^ ^ and W',W" corresponding walks in Ga,b constructed as in 
Proposition \5. li The projection 7r(C) of C is defined as follows: 

i. ifC is direct then 7r(C) = {W , W"}; 

ii. if C is twisted then tt{C) = W'U W" . 

Example 5.6. We return to the DSBP'^ and DSR graphs in Example \3. 2\ to illustrate the projections of 
direct and twisted cycles. The cycle 

(13,2^23,3^24,1^14,3^13) 
(shaded in the DSR^^^ graph on the left) is direct, and projects to the two closed walks 

W = (5i, i?2, 52, 5i), W" = (53, i?3, 54, i?3, 53) 
(shaded in the DSR graph on the right): 



I 

li) 



f 



31 



2^ 



23 



/ — 



i?2 



/ g 
S3 — i?3 — 54 



5i 



On the other hand, the cycle 

(13,1^23,2^12,2^13) 
(shaded in the DSR^'^^ graph on the left) is twisted, and projects to the closed walk 

W U W" = (5i, 52, i?2, 53, -R2, 5i) 
(shaded in the DSR graph on the right) 

.21,— "--,,0 9 
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52 
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S3 i?3 54 



Remark 5.7. Let C be a cycle in G[^^. 

L IfG IS twisted then \tt{G)\ = \C\', whereas if G is direct andTT{G) = {W',W"} then \C\ = \W'\ 



\W" 
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2. In the case that C is direct, various scenarios are possible: one of W or W" may be empty; W' and 
W" may be disjoint; or they may share vertices or edges; or they might even coincide. To see how one of 
W or W" may be empty, consider the following example: 



(12)--- 

An I 



A31 



(23) 



^i^i?i^ 



An 



c 



2^ 



W : 



W" 



25 —^—(24) 



i?4 



^5 



-i?3 



i?2 



- S'4 



Note that in this case, the index 2 is common to all S-vertices in C , and, disregarding signs of edges, the 
edge-sequence of W" is precisely that of C. In situations such as this, we make an abuse of terminology and 
write 7r(C) — W" rather than 7r(C) — {0,W"}. To see how W and W" may coincide consider the example 
illustrated below. For simplicity only S-vertices are depicted: 



C 




W',W": Si 



Identifying whether a cycle is direct or twisted. If ij and i'j' are two consecutive S-vertices in 
a cycle C in G\^ g, then {i,j} n {i',j'} must have exactly one element, say i. We call the pair {ij,ij') an 
inversion if {i — j){i — j') < 0. Whether a cycle is direct or twisted can be computed by counting inversions: 

[21 

Proposition 5.2. A cycle C of Gj^ g is direct if it has an even number of inversions, and twisted otherwise. 

Proof. Let C be a direct cycle with the sequence of S-vertices [ij — ai6i, a2b2, . . . , arbr = ij). Assume that 
Ofe < bk for all k. Recall that, by the construction in Proposition 15. ![ one element of {ak,bk) is assigned to 
W and the other is assigned to W" . Define md{akbk) as 



md{akbk) 



1, if Ok is assigned to W 

2, if Ok is assigned to W" . 



Note that uid{akbk) ^ ind(afc+i&fe+i) if and only if the pair (akbk, Ok+ibk+i) is an inversion. By definition, 
C is direct if and only if W' is a closed walk, i.e., if hid{arbr) = ind(ai5i), which is equivalent to saying that 
the number of switches between 1 and 2 in the sequence 



ind(ai6i), ind(a262), • ■ • , ind(ar6r) 



IS even. 



□ 



Example 5.8. The parity of the two cycles in Examvle \5.6\ can now be seen via the number of their inversions. 
The first cycle contains no inversions and is direct, whereas the second one contains one inversion (23, 12) 
and is therefore twisted. 

[21 

The parity of a cycle in Gj^ g relates to the parity of its projection in Ga.Bi as shown in the following 
result. 



22 

[21 

Proposition 5.3. Let C denote a cycle in Gj^ g. 

(i) IfC IS direct and n{C) = {W',W"} then P{C) = P{W')P{W"). 
(it) IfC is twisted then P{C) = -P{tt{C)). 

Proof. Let C = (ei, . . . , ca?) such that the initial vertex of ei is an S- vertex, and let {ii, . . . , im}, {ji, ■ • ■ , Jn-m} 
be the partition of {1, . . . obtained as in Proposition 15. II In view of Remark 15.41 we have 

M 
2 

(15) P{W') = (-1)^ n [sign(^(e.2._J)sign(^(e.,,_,+i))] . 

Equation ([14]) implies that, for any rG{l,...,^}, 

sign(7r(e,2^_J)sign(7r(ej2,__,+i)) = -sign(ej2^_ Jsign(e,2^_i+i) 

if the two S- vertices in the sequence e^j^ j , e^j^^j+i form an inversion, and 

sign(7r(e,2,_i))sign(7r(ej2,_i+i)) = sign(e,2^_ Jsign(ej2,_i+i) 

otherwise. Denoting by inv(VF') the number of pairs from (e^j , ei^+i), (e^g , 6^3+1), (eijj,^_^ , e^jj) whose 
S-vertices form an inversion, equation (jl5p becomes 

M 

(16) P{W') ^ (-l)f (-l)'-(^') nsign(e,,,_Jsign(e.,,_,+i). 

r=l 

Multiplying with the corresponding equation for W" yields 

N 

P{W')P{W") = (-1)*(-1)'"^('^')(-1)^(-1)'"^(^") ]^sign(e^) = (-1)'"^(<^)P(C) 

where inv(C) = inv(VK') +inv(VK") is the number of inversions in the S- vertex sequence of C. It follows from 
Proposition 15.21 that if C is direct then inv(C) is even and part (i) of the conclusion follows. If C is twisted 
then inv(C) is odd, P{tt{C)) = P(W' U W") = P{W')P(W") and (ii) follows. □ 

Having dealt in Proposition 15.31 with the behaviour of parity under projection we next ask what can be 

[21 

inferred about s-cycles in Gj^^ from their projections. 

[21 

Proposition 5.4. Let C be a cycle ofGj^g. 

(i) If C is twisted then tt{C) is an s-walk if and only if G is an s-cycle. 

(ii) If G is direct, tt{G) = {W ,W"}, and both W' and W" are s-walks, then G is an s-cycle. 

Proof. We may assume that G is directed; otherwise, if all its edges are unoriented, we simply choose one of 
the two possible orientations. By definition, G is an s-cycle if 

(17) n n '(e)- 

eG{S-to-R edges of C} eG{R-to-S edges of C} 

[21 

Since an S-to-R edge of G'^^'g projects to an S-to-R edge of Ga,b with the same label, one has 

n '(e) = n '(e) 

ee{S-to-R edges of C} ee{S-to-R edges of Tr(C)} 

in case (i), and 

n '(e) = f n '(e 

ee{S-to-R edges of C} \ee{S-to-R edges of W'} 
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in case (ii). The similar equalities hold for R-to-S edges and the conclusion is now immediate from (fTTl) . 

□ 

In conjunction with Lemma 12.191 Proposition 15 .41 has an immediate notable consequence: 

[21 

Corollary 5.9. If all cycles ofGA,B o,fe s-cycles, then the same is true for G]^ g. In particular, if Ga,b is 
acyclic, then every cycle of G^^ g is an s-cycle. 

Remark 5.10. The condition that all cycles in Ga.b are s-cycles is often fulfilled in applications to CRNs. 
Corollary 1 5. 9\ shows that when this is the case, verifying that G^^g is odd* boils down to ruling out odd 
intersections of e-cycles. 

5.3. Liftings of DSR cycles. As explained in Remark l5.3[ the projection map tt is not injective; it is easy 

[21 

to see how different cycles in Gj^ g may project to the same closed walk Ga,b- If is a cycle of Ga,b then 

any cycle G of G^^g that projects to W will be called a lifting of W. By way of Definition 15.51 and Remark 
[5Jl either G is twisted or, if 7r(C) = {W, W"} then W^' = or W" = 0. In the first case C wih be called 
an internal lifting of W and in the latter G is called an external lifting of W. The next proposition sheds 
light on this terminology. 

Proposition 5.5. LetW be a cycle inGA,B, o-nd let (Si-^, Rj-^^ . . . ^ Si^^ , Rj^^) denote its sequence of vertices. 

[2] 

(i) A cycle G in G^s '"^ '^'^ external lifting of W if and only if its vertex sequence is 

C= {iip,jf,...,iNp,jlf,iip) 

for some p ^ {ii, . . . , ijv}- 

(ii) If G is an internal lifting of W then the S-vertices of G are a subset of {ii, . . . , j^v} x • ■ • , iN}- 

Proof, (i) In one direction we can easily check that a cycle of the form G is an external lifting of W. On 
the other hand, suppose an external lifting of W gives rise to a cycle whose vertex sequence is not of this 
form, implying that G contains three consecutive S-vertices {ij,ij' ,i' j')^ with no common element (namely 
i,i',j,j' are all distinct). It follows from the proof of Proposition [5lT] that each of W and W" must include 
a vertex from {Si, Sj} and a vertex from {Si', Sj'}, which means that each of W and W" contains at least 
two S-vcrtices and therefore at least two edges. This contradicts the fact that one of them must be empty 
(ii) In this case, since G is twisted, neither of W or W" is empty. If ij is an S-vertex of G then, by 
construction. Si is an S-vertex of W and Sj is an S-vertex of W" . The conclusion follows. □ 

Remark 5.11. It is clear from the definitions and Provosition 1 5. 5\ that if a cycle G in G^ g is an external 
lifting of a cycle G' in Ga,b, then P{G) — P{G'), while if G is an internal lifting of G' , then P{G) — —P{G'). 
In other words internal liftings change parity, while external liftings preserve parity of cycles in Ga.b- 

Remark 5.12. A cycle G in Gj^ g of length 4 must be of the form {ap, cP , dp, , ap), with a ^ d and c ^ e. 

Thus, by observation 7r(C) is {Sa, Rc, Sd, Re, Sa), namely a cycle of length 4 in Ga,b. As G is clearly direct, 

it has the same parity as 7t{C). Thus by Provosition \5 . 5\ cycles of length 4 in G[^^ arise precisely as external 

liftings of cycles of length 4 in Ga,b- 

^ [2] 
Example 14.51 revisited. We may now argue that G^ j in this example is odd directly from Ga.i- Consider 

a cycle E of Gj^j, and let W' and W" denote the two walks in Ga.i that form its projection. Since there 

are (2) = 6 S-vertices in gJ^'j, the length of E can not exceed 12. Notice that Ga.i has two cycles: an even 
Hamiltonian cycle G, and an odd ^-cycle D. In view of Remark \5.4^ 1), neither W' nor W" can traverse 
back and forward consecutively an unoriented edge of Ga,i, o,nd t:{E) must therefore consist of copies of G 
and D. In fact, there is at most one copy of each of G and D in this union: two copies of G would imply the 
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contradiction lcngth(i?) > 16, while two copies of D would imply that E contains both liftings of D allowed 
by Remark \5. 12\ (where p is 3 and 4 respectively), so E cannot be a cycle itself. Now suppose that W' U W" 
has the same edges as C U D (counted with multiplicity). Then the directed path {S2, Si) is traversed 
twice, and therefore both directed paths Pi = (23, 2^, 13) and P2 = (24, 2"^, 14) mu.st exist in E. Further, 
{S4, i?4, 5*3) must arise as the projection of one of the directed paths P3 = (14, 4^, 13) or P4 = (24, 4^, 23). 
However, clearly P3 and Pi cannot both exist in E and similarly P4 and P2 cannot both exist in E, a 
contradiction. We conclude that E is either a lifting of C or a lifting of D. Since all liftings of e-cycle D 
are internal and all liftings of o-cycle C are external, it follows, by Remark \5.11[ that gJ^'j is odd. 



Remark 5.13. Although this could not occur in the previous example, in general the same cycle in Ga,b 
may lift to multiple cycles of different parities in g||'^, by lifting both externally and internally. This is 
illustrated in the following example: 



S2 



W : Ri 



S. 



Ri Si 
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i?2 Si 



A21 



Ail 



Ci : ^ 

A54 I 



14 

I A12 A24 
C2 : 

A52 A54 



, A21 
23 --- 



43 



I Ai 



{15}- 



W is an e-cycle, as is its external lifting C2 (note that it has four inversions). On the other hand, the 
internal lifting C'l is an o-cycle with one inversion. 

6. Further applications and limitations of the DSR^^l graph 

We present some examples which iUustrate use of the theory developed so far. In particular we show how 

it is possible in some cases, via the theory in Section [5l to make claims about (A_B)[^1 by examination of the 

[21 

DSR graph Ga.b without explicit construction of Gj^g. 



6.1. When the DSR graph has three S-vertices. If A G M^x, 

particularly nice structure, illustrated generically in the figure below: 



and B e Qq{A), then gJ^'^ ^as a 
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Depending on the particular entries of A and B, some edges may be absent, some edges may become directed, 

[21 

some edges may change sign, and some labels may become oo. Nonetheless, the generic structure of G^^ g is 
always the same and in particular, the degree of any R-vertex is at most 2. This fact implies that no odd 
intersection of cycles is possible (see Remark l2.16|) and checking that cj^'g is odd* amounts to showing that 
all its e-cycles are s-cycles. This has a pleasant translation into a property of the DSR graph, as follows: 

Theorem 6.1. Let A £ E^xm ^ ^ Qo{A*). Then, if Ga,b satisfies 

i. all e-cycles of length 4 are s-cycles; and 

ii. all o-cycles of length 6 are s-cycles, 
then G^"^ g is odd* , and {AB)^'^^ is a Po-matrix. 

Proof. We need to show that all e-cycles of gJ^'^ are s-cycles. As explained in Remark 15.121 all 4-cycles of 

[21 

B external liftings of 4-cycles of Ga,b, whose parity (and edge-labels) they inherit. It follows from 

(i) that all e-cycles of length 4 in G^^g are s-cycles. On the other hand, any 6-cycle of G^^g contains all 

S-vertices 12,23,13 of gJ^'^ and has therefore exactly one inversion. It follows from Proposition 15.41 that 
such cycle is an s-cycle if and only if tt{G) is a closed s-walk. 

We may distinguish three types of cycles, according to the three possible structures of their projections, 
illustrated below. Here vertices depicted as distinct are distinct, and the projected closed walks are obtained 
by following the arrows. Two arrows between a pair of vertices means that the corresponding (undirected) 
edge is traversed twice: 




Since the parity of a twisted cycle changes under projection, we need to show that odd walks in Ga.b of 
types (a), (b) and (c) are closed s- walks. Type (a) consists of three repeated edges and is structurally a 
closed s-walk (see Corollary I2.20p . Type (b) is composed of a repeated edge and a 4-cycle. If the walk is 
odd then the 4-cycle must be even and thus an s-cycle by hypothesis (i); it follows that the walk itself is a 
closed s-walk. Finally, type (c) is covered by hypothesis (ii). □ 

6.2. When the DSR graph is steady and either S- or R-vertices have degree no more than 2. 

The following generalisation of Proposition 14.21 is immediate: 

Theorem 6.2. Let A G ]R"X'" be such that Ga.a* is steady. 

(i) If A has no more than two nonzero entries in each column, then for any B £ Qa{A*'), (Ai?)[^l is a 
Po-matrix and so the nonreal spectrum of AB does not intersect C_. 

(ii) If A has no more than two nonzero entries in each row then the nonreal spectrum of AB does not 
intersect C_ . 

Proof. By CoroUarv 15.91 all cycles of G^^ are s-cycles, and the proof proceeds exactly as that of Proposi- 
tion 1121 ' □ 
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6.3. When the DSR graph is acyclic. 

Theorem 6.3. Let C be any n x m matrix such that Gc,c^ acyclic. Then, given any A G Qo(C) and 
B G Qo(C*)> is positive semistable. 

Proof. If the result holds for A G Q{C),B G Q(C*), then it follows for A G Qa{C),B G Qo(C*) since the 
set of positive semistable matrices is a closed subset of R"^". So assume that A G Q(C), B G Q(C*). Since 
is acyclic, it is odd, and so by Theorem 12 . 1 71 AB is a Po-matrix. Note further that since all edges in 
Gc,c* ^-re undirected, Gc,c*i being acyclic, must in fact be a forest. 

Now consider L for arbitrary A G Q(C), which clearly lies in Q(L ). Since G^.^* is acyclic, by 
Corollary 15.91 every cycle in G||'^t — Gj-a fj^^y is an s-cycle. By Lemma |2.21[ (AB)^'^^ = L^^L^ is a 
Po-matrix. Since both AB and (AB)!^! are Pg-niatrices, the conclusion now follows from Lemma □ 

Remark 6.4. Lemma \4.1\ is clearly a special case of Theorem \6.3[ If C is a column vector then Gq.c* is in 
fact a star. 

[2] 

Remark 6.5. // Gc.c* is a tree, this does not imply that is odd*. While most of the examples 

presented have involved situations where the DSR^^^ graph is odd or odd*. Theorem \ 6.3 illustrates that the 
D^i?!^] graph may provide useful information even in the case where it fails to be odd*. 

6.4. Limitations of the theory developed so far. There are a number of ways in which naive application 
of the theory described here can fail to provide information, even when a set of matrices structurally forbids 
the passage of a pair of nonreal eigenvalues through the imaginary axis. We illustrate with a couple of 
examples. Following up on Remark 16.51 the next example illustrates that even when gJ^'^ is not odd* and 
in fact (AP)[2l is not necessarily a Pg-matrix, it may still be possible to infer the absence of Hopf bifurcation 
from the DSR^^l graph: 

Example 6.6. Consider the matrices 

/ a 
A= b 
\ c 

[21 

with a, 6, c, d, e, f,g,h> 0. Gj^ g is shown: 




Labels have been omitted, but it is easy to check that Gy^ ^ fails to be odd* as a consequence of the cycle 
C — (13, 2"^, 23, 1"^, 13) which is even but, involving oo labels, not an s-cycle. Indeed (AB)^^^ is not necessarily 



/ ~d \ 
B=[ -e / , 
\ 9 h J 
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a Po-matrix for all values of the entries. However computation reveals that {AB)^'^^ is nonsingular and 
applying Lemma \8.5\ Spec (AB)\{0} does not intersect the imaginary axis, forbidding Hopf bifurcation. 



In fact, more careful application of the theory developed in |4j allows us to infer from G 



[2] 

A,B' 



and without 



direct computation of det{{AB)^^^), that det((^i?)f^l ) > 0; the argument is only sketched here. There are 
exactly two cycles in cj^'^, an odd cycle D = (12, 2^, 13, 1'^, 23, 3^, 12) and the even cycle C . However C does 
not contribute any negative terms to the determinant: this follows from the fact that no undirected edges are 
incident into S-vertex 12 and consequently there is no set of three R-vertices matched with the three S-vertices 



both G 



[2] 
A,0 



d clf'g in such a way that the union of these matchings includes C. See Remark ] 2. 18[ see 



also the discussion on "good" cycles [25| (resp. "central" cycles in and Pfaffian orientations. Since 
G is the only even cycle in G^^ g and contributes nothing to the determinant, det{{AB)^'^^) > 0. That 
det((Ai?)f^l ) ^ follows from the existence of D, which traverses each S-vertex and thus contributes a 
nonzero term to dct((v4B)[2l). 

Remark 6.7. Remark \2. 7| does not apply to Example \6.6\ since (AB)^'^^ has not been proved to be a Po-matrix. 
The next example demonstrates that (AB)^'^^ may be a Pg-matrix for all B S Qo{A*) although this does 

[21 . 

not follow from examination of G^ using any of the theory presented so far. It seems likely that this 

[21 

information can indeed be inferred from GJ.t, but that new techniques will be needed. 



Example 6.8. Consider the matrix 



A 



V 



\ 




1 

^1 / 



Symbolic computation confirms that (AB)^'^^ is a Po-matrix for arbitrary B G Qo(^*) • Since Ga.a* is odd. 



AB is also a Po-matrix, and so, by Lemma 

.[2] 



AB is in fact positive semistable. However, neither G 



A.A* 

[21 [21 

nor G it A is odd!' . The picture on the left shows two e-cycles of Gj^ j^t with odd intersection. On the right 



[21 

A* A depicted, where A denotes the submatrix of A obtained by removing its first and last rows. The 



-cycles (13,21,12,32,23,13,13) and (13, 3\ 14, l'^, 24, 3^, 23, 1^, 13) have odd intersection, and G?l ^ fails 



to be odd*. By Lemma \5.1[ Cj^l ^ 



is not odd* either. 



--,14, 



® 



43 




31 



l3 



12)-.., 2^ 
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Notice that although Gj^^t is steady, Lemma \2.21\ cannot be applied since G^l ^ is not in general steady 



An interesting feature of this example is that some products of the form contain 
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negative terms, but all such negative terms are cancelled by positive terms arising from other products. This 
may be understood as follows: given an edge-set E in Ga,b, '^^^{E) may include several matchings in G^^^, 
corresponding to terms of equal magnitude but possibly different sign in minors of [AB)^'^h the net result 
depends on the number of terms of each sign. 

The last two examples demonstrate that DSR^^l graphs have more to offer than is explored in this paper. 
In particular a more complete analysis of projection and lifting may provide new techniques to rule out Hopf 
bifurcation in families of matrices. 

We conclude by remarking that criteria such as "(A_B)[^1 is a Po-matrix" or even "det((Ai?)[^l ) > 0" which 
forbid the passage of a real eigenvalue of {AB)^^^ through zero are sufficient but not necessary to preclude 
Hopf bifurcation. Indeed a real eigenvalue of (AB) [■^l may pass through zero as a consequence of migration 
of eigenvalues of AB on the real axis without any possibility of Hopf bifurcation. This limitation may be 
circumvented by resorting to conditions that distinguish imaginary pairs of eigenvalues of AB from pairs 
of real eigenvalues summing to zero, for example based on subresultant criteria [101 . Natural extensions of 
the work here involve exploiting such criteria, and also examining other graphs related to AB whose cycle 
structure contains information about Hopf bifurcation. 
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